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A theory of the macroturbulent instability in the system containing vortices of opposite directions 
(vortices and antivortices) in hard superconductors is proposed. The origin of the instability is 
connected with the anisotropy of the current capability in the sample plane. The anisotropy results 
in the appearance of tangential discontinuity of the hydrodynamic velocity of vortex and antivortex 
motion near the front of magnetization reversal. As is known from the classical hydrodynamics 
of viscous fluids, this leads to the turbulization of flow. The examination is performed on the 
basis of the anisotropic power-law current-voltage characteristics. The dispersion equation for the 
dependence of the instability increment on the wave number of perturbation is obtained, solved, 
and analyzed analytically and numerically. It is shown that the instability can be observed even at 
relatively weak anisotropy. 
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I. INTRODUCTION 

Magnetic flux dynamics in type-II superconductors has 
been extensively studied since the end of the 50-s, start- 
ing with the pioneering work by A. A. Abrikosov. Pri- 
mary attention was given to hard superconductors, whose 
magnetic characteristics are defined by the presence of 
vortex pinning centers. The main features characterizing 
the nonuniform penetration of magnetic flux into such 
systems were revealed and studied; various theoretical 
models of electrodynamic processes in superconductors 
were suggested. An avalanche of new researches in this 
field was triggered by the discovery of high-T c supercon- 
ductivity (HTS), which brought into play thermal fluc- 
tuations and the strong anisotropy of superconductors. 
Different types of the phase transitions (melting, glass 
state transitions, etc. ) in the flux line lattice (FLL) 
have been discovered, studied, and explained. Many of 
the newly obtained results are described in comprehen- 
sive review papers^. 

The use of high-resolution magneto-optical (MO) tech- 
nique enabled an in-depth study of the dynamics of mag- 
netic flux in superconductors. One of the most impor- 
tant features revealed by means of this method is the 
fractal structures of thermally activated flow of the mag- 
netic flux. Such fractal structures arise usually due to 
the development of the characteristic instabilities such as 
macroturbulence in 1-2-3 systems 3 ' 4 ' 5 and dendrite insta- 



bility in magnesium deboride 6 . However, these dynamic 
instabilities in FLL have not been so far properly inves- 
tigated. 

Perhaps, the macroturbulence is the most dramatic 
phenomena observed in the dynamics of the magnetic flux 
in HTS. It appears like a turbulization of the FLL motion 
in a sample near the front of magnetization reversal that 
separates regions of the vortices of opposite directions 
(vortices and antivortices). Note that the macroturbu- 
lence was only revealed in single-crystal samples of the 
1-2-3 system. Its essence is as follows. When a magnetic 
flux is trapped in a superconductor and a moderate field 
of a reverse direction is subsequently applied, a bound- 
ary of zero flux density will separate the regions contain- 
ing vortices and antivortices. For definiteness, we apply 
the term "antivortices" to the vortices whose direction 
coincides with that of the external magnetic field, and 
the term "vortices" , to the vortices which were originally 
present in the sample, prior to switching on the magnetic 
field of negative sign. At some range of magnetic fields 
and temperatures, this flux-antiflux distribution becomes 
unstable. A disordered motion of magnetic flux arises 
at the front of magnetization reversal, which resembles 
a turbulent fluid flow. This process rapidly develops in 
time and is accompanied by the emergence of channels via 
which the antivortices penetrate into the region occupied 
by the vortices. In other words, the front of magnetiza- 
tion reversal takes on a finger-like shape. The annihila- 
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tion of vortices and antivortices occurs at the front, and 
the process of macroturbulence is soon terminated after 
a complete disappearance of the vortices. This pattern 
of penetration of the magnetic flux differs qualitatively 
from the steady-state slow motion of the front of magne- 
tization upon initial switching on of the magnetic field, 
when vortices of only one direction are present in the 
sample. Note that the characteristic times of instability 
development amount to several seconds and more, and 
emerging spatial structures are macroscopic, i. e., they 
contain a large number of single vortices. 

The described instability cannot be understood within 
the framework of the critical state model or conventional 
models for the thermoactivated flux relaxation. At the 
same time, this phenomenon is a close analogue to the 
turbulence in the usual hydrodynamics. So, the study 
of the phenomenon is of general physical interest. More- 
over, this investigation has evident application aspects, 
since the macroturbulence can affect the ac losses, mag- 
netic noise, and relaxation phenomena in superconduct- 
ing devices. 

An attempt to explain this remarkable behavior of the 
flux-antiflux interface was undertaken in paper by Bass et 
al. — , where the instability was attributed to a thermal 
wave generated by the heat release due to the vortex- 
antivortex annihilation. Unfortunately, this mechanism 
can hardly be relevant. Indeed, the annihilation energy 
is small and a corresponding temperature rise is negligi- 
bleS& 

Another physical pattern for the emergence of macro- 
turbulence was discussed by Vlasko-Vlasov et al—. They 
draw attention to the fact that the process of annihila- 
tion of a vortex-antivortex pair may be accompanied by 
the formation of spatial domains free of vortices (the so- 
called Meissner holes) . It was assumed that the presence 
of such domains may cause instability due to high cur- 
rents which have to flow in the vicinity of such a Meissner 
hole. Yet, the authors of Ref.^did not describe the pos- 
sible physical instability mechanism. 

An explanation of the macroturbulence should focus 
on the experimental fact that the instability was reported 
for YBa2Cu307_5 and other 1-2-3 single crystals only, 
which are characterized by the anisotropy in the basal ab 
plane. This anisotropy can be related to a specific crys- 
tallographic structure of these HTS and, in particular, 
to twin boundaries^. The electromagnetic instabilities 
of the critical and resistive states in anisotropic super- 
conductors were considered by Gurevichi2*i£. However, 
these results cannot be directly applied to the explana- 
tion of the turbulent instability. 

An alternative approach to understanding the mecha- 
nism of the macroturbulence was elaborated in Refs. 00 
taking into account of the specific features of the flux mo- 
tion in the anisotropic superconductors. The anisotropy 
gives rise to the motion of the flux lines at some angle 
with respect to the Lorentz force direction. For example, 
in the presence of twin boundaries, vortices and antivor- 
tices move preferably along these guiding boundaries (the 



so-called guiding effeclii^ii^) . In general, the angle be- 
tween the twins and the crystal grains is around 45°. As 
a result, the flux lines move at some angle with respect 
to the magnetization reversal front. In our opinion, it is 
just this circumstance that is of paramount importance 
to ascertain the nature of macroturbulence. The vor- 
tices and antivortices are forced to move towards each 
other in such a way that the tangential component of 
their velocity becomes discontinuous at the flux-antiflux 
interface. According to a classical paper of Hclmholtz, 
a stationary hydrodynamical flow can be unstable and 
turbulent under such conditions. It was shownS that a 
purely hydrodynamic approach which takes into account 
the anisotropic viscosity of the flux flow can provide the 
basis to gain an appropriate insight into the nature of 
the macroturbulence. In particular, the macroturbulence 
is usually observed within a rather narrow temperature 
window in the vicinity of 40-50 K—. As was demon- 
strated in Rcf. |3, the temperature window is consider- 
ably wider for heavily twinned samples in which a more 
pronounced anisotropy can be expected. 

The model developed in Refs. 00 is based on the anal- 
ysis of the viscous magnetic flux flow in anisotropic su- 
perconductors. In the framework of this model the de- 
pendence of the viscosity coefficient of the flux line lattice 
on the flux velocity was neglected. This approximation 
corresponds evidently to the description of the supercon- 
ductor in terms of the linear current- voltage (I—V) char- 
acteristics. This idealized model allows one to describe 
qualitatively the effect in question. However, it predicts 
that the macroturbulence is observable for unrealistically 
high values of the anisotropy parameter. Note that the 
regime of the viscous magnetic flux flow is realized at ex- 
tremely high current densities that cannot be achieved 
under experimental conditions used in the MO measure- 
ment. This is a reason for using a more general model for 
the flux flow in the present paper. Specifically, we take 
into account the dependence of the viscosity coefficient 
on the flux flow rate. In other words, we describe the 
electric properties of a superconductor by a more realis- 
tic non-linear I—V characteristics. We have theoretically 
studied the instability of the magnetization reversal front 
under the assumption that the current-voltage character- 
istic is a power-law function with an exponent m > 1. We 
show that, even at a comparatively weak anisotropy of 
the current-voltage characteristics, the flow of a system 
of vortices and antivortices in a superconductor can be 
unstable. A preliminary results of the present study was 
published in the letter—. 



II. MAIN EQUATIONS 

The macroturbulence is usually observed in plate-like 
single crystals placed in a transverse magnetic field. The 
MO image provides information about the distribution of 
the normal component of magnetic induction. The pene- 
tration of an electromagnetic field into a superconductor 
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in such a geometry is of interest per se and was studied 
by numerous researchers (see, for example, Rcfs. 18 - 19 ). 
However, the turbulent behavior of magnetic flux is not 
a geometric effect. It was observed in thin plates and in 
single crystals with a low demagnetizing factor as well. 
Frello et al£ reported the MO visualization of developed 
macroturbulence of vortex matter in an Nd-123 crystal 
3.1 x 2.5 x 1.3 mm 3 in size. In the magnetic field, this sam- 
ple was divided into three magnetically unbound regions 
(each having smaller transverse dimensions) in which the 
turbulence was developing independently. Moreover, it 
should also be kept in mind that the instability is of- 
ten observed under conditions of the full penetration of 
magnetic flux into the sample when the difference in the 
distribution of induction in the cases of longitudinal and 
transverse geometry turns out to be not too significants^. 
Therefore, for a qualitative description we can restrict 
ourselves by the simplest longitudinal geometry. So, con- 
sider an infinite superconducting plate of thickness 2d in 
the external magnetic field H directed parallel to the 
sample surface along the z-axis. The cc-axis is perpendic- 
ular to the plate and the x-axis origin, x — 0, is placed in 
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median plane x = of the plate. 

In real crystals, the vortices and antivortices are pinned 
by defects. The vortex motion occurs due to the action of 
the Lorentz force and (at non-zero temperature T) due to 
the thermoactivated flux creep. In the framework of the 
macroscopic description of the flux flow we can take this 
effect into consideration by making use of the non- linear 
I —V characteristics 1 . There exists an additional reason 
for the flux line motion in the situation under study. The 
vortices and antivortices can annihilate near the plane of 
zero magnetic induction. This phenomenon gives rise to 
diminishing the vortices number in the bulk and to en- 
tering new antivortices through the sample surface. As 
a result, the line of zero induction moves towards the 
sample center with time. We will use the phenomenolog- 
ical (hydrodynamic-like) equations to describe the flux 
line motion. This approach is valid if all spatial scales 
of the problem are much larger than the FLL constant 
df. In particular, the characteristic spatial scale of the 
macroturbulence l c should be large, l c ^> df. 

Let us denote the density of the vortices and antivor- 
tices as Ni(x,y) and N 2 {x,y), respectively. The relation 
between vortex densities N a (x,y) and the magnetic in- 
duction B(x,y) is evident, 



N a (x,y) = s a B(x,y)/$o, 



a = 1,2 



(1) 



where si = 1, s 2 = — 1 and <!>o is the magnetic flux 
quantum. ^From the symmetry of the problem, it is suf- 
ficient to consider only the region < x < d. Fig. [21 
shows schematically the spatial distributions of Ni(x) 
and N2(x). The vortex and antivortex densities should 
evidently satisfy the continuity equation, 



dNg 
dt 



+ div(iV Q U Q ) = 0, 



(2) 



FIG. 1: Geometry of the problem. 



where V a denotes the hydrodynamic velocities of the vor- 
tices and antivortices. 

In the hydrodynamic approximation, there are two 
equivalent approaches to find the vortex and antivor- 
tex velocities. The usual dynamic approach, where V a 
is defined by the viscosity equation for FLL, which t akes 
account of the Lorentz force, was applied in Refs. I8l9l 
The alternative approach implies the use of the I — V 
characteristics, 



field H first increase up to some value much higher than 
the lower critical magnetic field H&. Then H is lowered 
through zero to a negative constant value whose modulus 
is likewise higher than H c \. Under such conditions, two 
kinds of vortices exist inside the sample. The first one 
with the magnetic flux directed along the positive z-axis 
has penetrated into the sample under the magnetic field 
rise [vortices). The vortices occupy the central part of 
the superconductor. The second group of vortices en- 
tering the sample after the external field has changed 
its sign is characterized by the flux directed oppositely 
[antivortices). The antivortices are located in two pe- 
ripheral regions arranged symmetrically relative to the 



J = J(E), 



(3) 



where J and E are the macroscopic current density and 
the electric field. We operate with the latter method in 
this paper. 

The velocities and electric field are interrelated by the 
usual equation, 



E = -V x B. 

c 



(4) 



Using Eq. Q), we can rewrite Eq. (0} in the form, 
N a s a § „ „ N a s a <$> 



E r 



-Vcty, E y 



V ax . (5) 
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FIG. 2: Flux distribution in one half of an infinite slab (|x| < 
d) containing trapped vortices of density Ni (x) in the central 
region \x\ < xo, and antivortices of density N2(x) penetrating 
from the outside. The other symbols are defined in the text. 



In principle, we can use the I —V characteristics in gen- 
eral form ©. However, in so doing the obtained results 
are rather cumbersome and their physical understand- 
ing is not evident. To avoid this inconvenience, we use 
here some particular power-law I — V characteristics for 
anisotropic superconductors. Let the main anisotropy 
axes be X and Y (e. g., the X-axis is directed along 
twins, while the Y-axis is across them) . We suppose that 
the components of the current density vector along these 
directions are the odd functions of the electric field and 
for positive E can be presented as 



e \E 



1/m 



l/m 



(6) 



where J c is the critical current density along the Y axis, 
to > 1 is the corresponding exponent, e < 1 is the 
anisotropy parameter for the critical current components 
and the value of J c is defined as the current density Jy at 
E = Eq (usually, Eq is accepted as 1 /zV/cm). At rn = 1 
the I — V characteristics © correspond to the viscous 
flux flow regime used for the macroturbulence analysis in 
Ref. In the limiting case m>l, these characteristics 
are in accordance with the critical state model proposed 
by ClemSl for the description of the electrodynamics of 
anisotropic superconductors. 

As was mentioned above, the anisotropy in ab plane 
arises in 1-2-3 single crystals mainly due to the exis- 
tence of twin boundaries. In a real experimental situ- 
ation, these boundaries make angles of about 45° with 
the axes x and y. Therefore, we assume the angle be- 
tween two coordinate systems, xy and XY, to be equal to 
45°. Now, using Eqs. JSJ, I© and the Maxwell equation 
VxB = 4irJ/c we can derive the equations interrelating 
the vortex density and velocity components, 



dN a dN a _ 4irV2J c 
dx dy c$o£ 



cEoV? 



(7) 



8N a 8N a 



dx 



dy 



(V xa + V ya ) 



7'. 



To solve the problem, we must formulate the bound- 
ary conditions at the sample surface, as well as at the 
interface between the domains occupied by the vortices 
and antivortices. We will first discuss the conditions at 
the sample boundaries. Ignoring the induction jump on 
the surface (this may be done in the case of fairly high 
values of H , H ^> H c i), one can derive 



N 2 (d) =N 2 (-d) =H/<Z> . 



(8) 



Since only the right-hand part < x < d of the sample 
is treated, we will replace the condition N2(~d) = H/$q 
by the requirement, 



14(0) = 0, 



(9) 



that immediately follows from the symmetry of the prob- 
lem. 

Now we turn to the boundary conditions at the moving 
vortex-antivortex interface. In general case, the position 
of the interface depends on time t and the y coordinate. 
The equation for this surface can be presented as x = 
xa(y,t). Then, one can define the interface velocity U as 
a vector normal to this surface, 



U x = 



dxn 



1 



dt l + (dx /dy) 2 



Uy = - 



dx a dxp/dy 
dt l + {dx Q /d y y 



(10) 



In general, the velocity U, just as xo(y,t), depends on 
time and the coordinate y, U = U(y, t). 

In the coordinate system which moves with the inter- 
face, the total flux of the vortex and antivortex densities 
through the interface should vanish due to the evident 
vortex conservation law. This means that the following 
boundary condition should be fulfilled at x = Xq, 



NxlVi-U) +N 2 [V 2 -U 



0, 



(11) 



where the subscript n denotes the vector components 
transverse to the interface. The components of the cor- 
responding normal unit vector v are 

1 



i + (dx Q /d y y 



dx /dy 



1 + {dxo/dyf 



(12) 



The second boundary condition should define the rate 
of the vortex-antivortex annihilation. It is obvious that 
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this rate goes to zero if the density of vortices or an- 
tivortices at the interface between them is zero. Then, 
following the conventional approach to describing such 
kinetic processes, we represent the rate of annihilation to 
be proportional to the product of the vortex and antivor- 
tex densities, 



Ni(Vi - U) n = RN ± N 2 . 



(13) 



A similar model for the annihilation process was used 
in Ref. 0. The finite size region where the annihilation 
occurs was supposed to exist in the sample. The annihi- 
lation rate was assumed to be proportional not only to 
the product N\ and N 2 but also to the relative velocity 
\Vi — V 2 \ of the vortices and antivortices. Contrary to 
Ref. 0, we believe that the region where vortices and an- 
tivortices coexist and annihilate is small enough and is 
much less as compared not only to the sample sizes but 
also to the characteristic spatial scale l c of the macro- 
turbulcncc. This assumption can be confirmed by MO 
images of the macroturbulent instability 3,4,5 . It means, 
in particular, that the annihilation constant R is not too 
small. Thus, the flux-antiflux boundary can be repre- 
sented by the surface x = Xo(y,t). 

Finally, we assume that the average magnetic induc- 
tion in the neighborhood of the interface is zero, i.e., 



N ± = N 2 



(14) 



at x = xo{y,t). One can readily demonstrate that this 
condition directly follows from Eq. J2J) and relation (|ll|l , 
if it was valid at the moment of emergence of antivortices 
into the sample. In our case, N% = N 2 = at the time 
moment when the decreasing external magnetic field as- 
sumed the value of H = 0. 



III. UNPERTURBED PROFILE OF THE 
VORTEX DISTRIBUTION 

The formulated set of equations and boundary con- 
ditions allow us, in principle, to analyze the behavior of 
the vortex-antivortex system. We assume that the unper- 
turbed flux-antiflux boundary is a straight line. A scheme 
of solving the problem consists in finding the unperturbed 
vortex distribution N a (x,t) and velocity V a (x,t). Using 
these base distributions, one should solve a linearized set 
of equations for fluctuations 5N a {x, y, t) and 5V a (x, y, t). 
This section is devoted to the analysis of the base profile. 
Below we take into account the dependence of the critical 
current density J c on the magnetic induction <f>aN a . For 
definiteness sake, we choose it in the simplest form, 



J c = A/N a . 



(15) 



To find the base profile of the vortex density it is nec- 
essary to solve a complex set of partial differential equa- 
tions © -0 along with the boundary conditions iJS}- 
(|14|l some of which are written at the moving interface. 



A peculiar complication consists in the fact that the ve- 
locity U (t) of the motion of the vortex-antivortex inter- 
face is unknown and should be found self-consistcntly. 
Unfortunately, the problem does not have simple (with 
U = const 7^ 0) automodel solutions because the in- 
terface moves non-uniformly under the conditions being 
studied. 

As the first approximation, let us calculate a station- 
ary base profile of another system where the velocity U(t) 
equals zero. A superconducting slab carrying the critical 
transport current is an example where such a distribution 
of the vortices and antivortices can be realized. Assum- 
ing dN a /dt = in Eq. (J2J, one can easily obtain the 
distributions N\(x) and N 2 (x), 



N a (x) = Nay/1 + s a C(x 



C = 



8ir\/2dA 
c$ A 2 



cE Q (l + e m ) 



x)/d, 

l/m 



(16) 



from Eqs. ©-0 and conditions l(TT )H fTi )l . Following 
the experiments^* 5 - which show that the vortex density 
Nq at the interface is small with respect to its value 
N 2 (d) = at the sample surface, we assume the 

constant C to be much higher than unity. Neglecting 
unity in the radicand in Eq. I|16(l and using condition |jSJ 
one can evaluate the vortex density Nq at the interface, 



N a = N a (x = x a ) ~ ( — 



2 (m+l)/2 $oi? 



1/2 



(17) 



Expression (|16J) is the solution of the stationary prob- 
lem. We are interested in the base profile corresponding 
to the moving interface. This motion leads to a distor- 
tion of the profile with respect to Eq. i|16|) . However, we 
assume the velocity U to be small with respect to the 
vortex velocity V a , U <C V a . As it can be readily shown, 
in this case the base profile differs but slightly from the 
stationary distribution (|16f) . 

It is convenient to perform the further analysis using 
the dimensionless variables, 



n a = N a /N , t = t/to, t = 
t = x/L, C = y/L, L = 

V2cE '' 



8irAE ' 

inVZA 1 
u = Uto/L. 



(18) 



The normalization to the time-dependent value = 
N a (xo(t)) is allowable here since we assume that the in- 
stability develops much faster than noticeable changes in 
N a (xo) and U(t) occur. 

Taking into account the smallness of the velocity u(t) 
of the flux-antiflux boundary, we can linearize the ex- 
pression for the base profile n a (£, t) with respect to u. 
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In so doing, we get the equations for the first and sec- 
ond derivatives, n a and n a , of the unperturbed vortex 
density near the interface, 



n a (x = x ) = -s a p ( 1 + s a 



n a (x = .T ) = -{n a f 



2u 



(l + e)m«) 



m-2 ' 



(19) 
(20) 



c = s 



P 



2r 
1 + e 



l/n 



These formula are used in the next section devoted to the 
analysis of the instability. 



2ik- 



(24) 



We assume the perturbation of the vortex density to be 
damped away from the interface £ = £o( r ) a t short dis- 
tances with respect to the sample thickness d. This allows 
us to replace the coordinate-dependent variables n Q (£) 
and fi a (£) in Eq. H24[) by their values n a (£ — £o) and 
n Q (£ = £q) at the interface (see Eq. jTSJ)). 

The solution of Eq. I|24|) defines the exponential behav- 
ior of the fluctuations /i and f2, 



IV. INSTABILITY IN THE 
VORTEX- ANTIVORTEX SYSTEM 



A(£-£o) = /i exp[ Pl (£-e )], 



Let the dimensionless vortex density be a sum of the 
unperturbed term n Q (£,r) and a fluctuation term, 

n a = h a + f a (£ - £o(T"))exp(Ar + ik(). (21) where 

The linearized boundary conditions should be written at 
the perturbed interface, 



C = ^o(C,r)=Co(T)+^exp(ifcC + Ar). (22) 
It follows directly from Eq. (f L L|t that 

5£ = (/i - h)/2p. (23) 



Substituting Eq. (|^T|) into Eq. and neglecting the 
terms proportional to e 2 , we derive the ordinary differen- 
tial equation with the coordinate-dependent coefficients 
for the fluctuation J a (£ — (,o(t)), 



ik(l - 2e) + ^ - u 

n a mha n r r 



2A 



~'(m— ~ ' (m — 2) ~ (m+1) 



2u 



2uik(m — 1) 

~'(m— l) ~ m 



/ 2 (£-£o) = / 2 exp[p 2 (£-&)], 



Px >2 = ±p - ik + pu/2r ± fi 1)2 , 



(25) 



(26) 



n x , a = 



2p 



ifcitp Xp 



if. H/' /v/' ,9 ikupim — l)u 
' -t-+ 4efc 2 =F 9 o 



± 



mr mr 



2-u/? 2 \p(m — l)u 



mr 



T 



1/2 



ReOi, 2 >0. (27) 



Here the upper and lower signs correspond to the coeffi- 
cients pi and P2 , respectively. 

Substituting Eqs. into the boundary condi- 

tions qilfl. iQ gives two linear algebraic homogeneous 
equations for /i and / 2 , 



+ ifc(l - 2e) + - iMl ~ 3c) - u - A 

P P 



• ~ ' (m — 1) 



(p 2 + ifc(l - 2e) + n 2 ) + 



i/cu(l — 3e) 



A 

u+ - 
P. 



= 0, 



(28) 



TO _'( m -l) 



(pi+iA(l - 2e) +%) + 



ife(l-2e)_> 
n-, 



4p 
2^~ 



A 

27 



2p 



ifc(l-2e) 
4^ 



u 
m 



2p 



0. 



(29) 
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Here we took into account that the normal unit vector to 
the perturbed interface is v = (1, — iM£(C,r)). 

Requiring that the determinant of set (j28(l , 129J1 should 
vanish and omitting the terms of the order of u 2 <C 1, 
one obtains the dispersion equation for the increment A 
at different values of the wave number k, 



A = - Aek z - iku— - 2p' 

p V mr 



Here fi is the root with Re ft > of the equation 
771 mm 

ik p(n 2 ^l + n pm ■ - 2 



(30) 



2 2 



4eAT ( fr + 0— + iku 

m 2m z r 



p 777 

, • (31) 



Below we perform the analysis of this equation in two 
qualitatively different limiting cases. 



A. Linear viscosity, m = 1 

The vortex motion in the case m = 1 can be de- 
scribed in terms of the magnetic flux viscose flow with 
the constant anisotropic viscosity coefficient. Indeed, the 
current- voltage characteristics (JHJ, (|15|l leads to the lin- 
ear tensor relationship between the vortex velocity and 
the Lorentz driving force, 
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-B x J, 
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where the viscosity tensor Tn~ is characterized by the 
principal values 
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At to = 1, Eq. 131(1 takes on the form, 

n 4 + 6rn 3 - 32r 3 fl - 2ikur(fl + 4r) - 
4ek 2 n(fl + 4r) = 0. 



(34) 



It has a simple solution with the following asymptotics 
at k — * oo, e — > 0, 



ft w (2fc|7t|r) 1/3 exp(-i7r/6). 



(35) 



The corresponding increment A increases unrestrictedly 
at k — > oo, 



A w (fc|7i|r/4) 2/3 - ifci 



(36) 



The existence of the solution with Re A > means 
that the instability of the base profiles near the flux- 
antiflux interface exists. The increment of the instability 
grows with an increase of the wave number if k < k m . 
Therefore, the perturbations with the characteristic spa- 
tial scale l c ~ l/k m along the y-axis is predominant over 
others. Thus, the flux-antiflux interface is disturbed in a 
turbulent-like manner— and we can suppose that the pre- 
sented mechanism of the macroturbulence could explain 
the experiment. 

As was shown in Ref. |8|, the instability exists for very 
small e only, 



e < 0.019 



U 



2RN 



< 1. 



(38) 



So stringent a requirement imposed on the parameter 
of anisotropy prompted us to consider the more general 
situation of arbitrary m. The case m ^> 1 is discussed in 
the next subsection. 



B. Power-law current-voltage characteristics, m 3> 1 

The parameter e = e m in Eq. I|31|) becomes negligible 
at 777 3> 1. Therefore, Eq. l|31|l transforms into 
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X 2 + X + 2 = 0, (39) 
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P 

Contrary to Eq. I|34|) . this equation contains the term 
iKX 2 (m — l)/m which plays an essential role at k > 1 
and changes radically the character of the solution: 

The correspondent increment A, 
A fts mrp I in + 2 j- 



2V2 



777, K>1, (41) 



is almost an imaginary quantity but its relatively small 
real part is positive at k < 8m 2 . This means that the 
instability develops as an oscillating process with a mag- 
nitude increasing in time. The wave number n of the 
instability is bounded from above due to the finiteness 
of the parameter l/m. Contrary to the case m — 1, the 
instability occurs regardless of the strength of the cur- 
rent anisotropy, practically, if the anisotropy parameter 
e < 1. 



Considering the finite value of e limits this increase to a 
certain maximum value of RcA m 8 , 



V. NUMERICAL ANALYSIS AND DISCUSSION 



On the basis of the developed theory, let us analyze the 
instability of the vortex-antivortex system. In general 
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FIG. 3: The numerical solution ReA(fv) of the dispersion 
equation for m = 20, (u/r) 2 — 0.02 at different values of the 
anisotropy parameters e: e — 0.5 (curve 1), e = 0.45 (2), 
e = s c = 0.43 (3), e = 0.4 (4), e = 0.38 (5), and e = 0.2 (6). 



case of arbitrary e and m this analysis requires that the 
dispersion equation (|31[1 be numerically solved. 

The dependence of the increment ReA on the wave 
number k for m = 20 and the different values of the 
anisotropy parameter e is shown in Fig. [3] For definite- 
ness, the value of the ratio (u/r) 2 is taken to be equal 
to 0.02. Remind that this ratio should be small in ac- 
cordance with the assumptions made above. The figure 
shows that the spectrum of perturbations strongly de- 
pends on the parameter of anisotropy. Instability is ob- 
served at e < e c ss 0.43. Indeed, the increment ReA 
becomes positive at some interval of the wave number 
k and reaches the maximum value at finite k. If the 
anisotropy parameter is higher than e C) the increment 
ReA is negative at any value of k. This implies that the 
vortex-antivortex system is stable with respect to small 
perturbations if the current anisotropy of a superconduc- 
tor is not small enough. Of course, the critical value e c of 
the anisotropy parameter depends strongly on the expo- 
nent m. The graphs in Fig. ^illustrate the change in the 
value of e c with an increase of m. The function s c (m) 
is not very sensitive to the value of the small parame- 
ter u/r in our theory. The graph £ c (m) represents the 
separatrix dividing the phase space (m, e) into two parts 
corresponding to the stable (left-hand part) and unstable 
(right-hand part) states of the vortex system. 

One of the most important results of this paper con- 
sists in the substantial weakening the requirement theo- 
retically imposed on the anisotropy parameter e for the 
observation of the instability in the vortex-antivortex sys- 
tem. The effective parameter of the anisotropy determin- 
ing the existence of the instability is e m in the consid- 
ered case of the nonlinear current-voltage characteristics. 
Since the exponent m can reach several tens in real super- 
conductors, the necessary condition for the instability is 
achieved for not too small values of e. This circumstance 
allows one to eliminate a seeming contradiction between 
our theory and experiments where the macroturbulence 




FIG. 4: The separatrix dividing the stability and instability 
regions in the phase plane (e,m). The solid line corresponds 
to (u/r) 2 — 0.02, the dashed line is obtained for (u/r) 2 = 
0.002. 



was reported to be observed even in detwinned YBCO 
single crystals. The MO image of the development of 
the macroturbulence in such a sample is shown in Fig. [3] 
The twinning structure in this sample is not observed 
in polarized light whereas the macroturbulence is clearly 
pronounced. However, a careful scan of the left top image 
in Fig. [5] where the initial magnetic flux distribution is 
shown, highlights the presence of some anisotropy of the 
current-carrying capability in the sample plain. The pen- 
etration depth in the horizontal direction is clearly seen 
to be about 1.2 times higher than in the vertical one. 
This means that the anisotropy exists and the macro- 
turbulent instability may appear. Such a situation with 
the detwinned samples is, perhaps, typical for the 1—2- 
3 single crystals. After the detwinning procedure, some 
traces of the twin structure remain inside a sample and 
the anisotropic distribution of the impurities takes place. 
The existence of the anisotropy of the electrical resistiv- 
ity in the ab plane of detwinned YBCO single crystals 
was reported in Ref. |2^ Irrespective of its nature, the 
anisotropy in the detwinned samples can cause the ob- 
served macroturbulent instability. 

Finally, we should note that the qualitative features of 
the studied instability are quite similar at any m > 1. 



VI. CONCLUSION 

The very interesting phenomenon of the macroturbu- 
lent instability has so far been observed only in super- 
conductors of 1-2-3 systems. Such systems are nor- 
mally characterized by a well-pronounced anisotropy of 
the current-carrying capability in the ab plane. This 
experimental fact provided the basis for the theoretical 
approach to explain the nature of the instability. Under 
the action of the Lorentz force, the vortices and antivor- 
tices move towards each other at some non-zero angle 
with respect to the front of magnetization reversal owing 
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FIG. 5: The MO image of the development of the macrotur- 
bulence in a detwinned YBCO single crystal. 

to the anisotropy. As a result, the vortex and antivortex 
hydrodynamic flow is characterized by the tangential dis- 
continuity near the front. Specifically this fact leads to 
the turbulization. The anisotropy is described in terms of 
the power-law anisotropic current-voltage characteristics. 
The analysis of the linearizes set of equations consisting 
of the continuity and Maxwell equations along with ap- 
propriate boundary conditions allowed us to derive the 
dispersion equation for the increment of the instability. 
As is shown, the instability exists in a wide range of the 
parameters of the problem. In particular, the instabil- 
ity can be observed at not very strong anisotropy if the 
current-voltage characteristics is steep. This conclusion 



agrees with the experiment. 

Unfortunately we cannot make a direct comparison 
between the theoretical results and the experimental 
data because the theory operates with the unknown 
phenomenological parameter R describing the vortex- 
antivortex annihilation rate. In order to express this 
parameter in terms of observable variables, the micro- 
scopic model of the annihilation process needs to be con- 
structed. Nevertheless, the proposed theory and the ob- 
tained phase diagram, showing the region where the in- 
stability occurs, can be useful for the qualitative analysis 
of the macroturbulence. For instance, the existence of 
the instability in a definite temperature region alone can 
be simply rationalized within our model. At low tem- 
peratures, the critical current density increases, and the 
characteristic spatial scale L, in Eq. (|18|) . decreases corre- 
spondingly and becomes comparable to or less than the 
twin-boundary spacing. As a result, the anisotropy is 
suppressed and the instability disappears. On the other 
hand, at temperatures close to T c the anisotropy is no 
longer effective due to the thermal activation of the vor- 
tices. It is a remarkable thing fully supporting our model, 
that in the present heavily twinned crystal the turbu- 
lence occurs at much lower temperatures than is found in 
previous studies of similar crystals with only little twin- 
ningvii2£. 
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